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ABSTRACT 

The  ultimate  goal  of  the  prediction  of  Sea  Surface  Temperature  (SST)  from  satellite  data  is  to  attain  an  accuracy  of 
0.3°K  or  better  when  compared  to  floating  or  drifting  buoys  located  around  the  globe.  Current  daytime  SST  algorithms 
are  able  to  routinely  achieve  an  accuracy  of  0.5 °K  for  satellite  zenith  angles  up  to  53°.  The  full  scan  swath  of  VIIRS 
(Visible  Infrared  Imaging  Radiometer  Suite)  results  in  satellite  zenith  angles  up  to  70°,  so  that  successful  retrieval  of 
SST  from  VIIRS  at  these  higher  angles  would  greatly  increase  global  coverage.  However,  the  accuracy  of  present  SST 
algorithms  steadily  degrades  to  nearly  0.7°K  as  the  satellite  zenith  angle  reaches  70°,  due  mostly  to  the  effects  of 
increased  atmospheric  path  length.  We  investigated  the  use  of  Tfield,  a  gap-free  first  guess  temperature  field  used  in 
NLSST,  as  a  separate  predictor  to  the  MCSST  algorithm  in  order  to  clearly  evaluate  its  effects.  Results  of  this  new 
algorithm,  TfieldSST,  showed  how  its  rms  error  is  heavily  dependent  on  the  aggressiveness  of  the  pre-filtering  of  buoy 
matchup  data  with  respect  to  Tfield.  It  also  illustrated  the  importance  of  fully  exploiting  the  a  priori  satellite-only 
information  contained  in  Tfield,  presently  tamed  in  the  NLSST  algorithm  due  to  the  fact  that  it  shows  up  as  a  multiplier  to 
another  predictor.  Preliminary  results  show  that  SST  retrievals  using  TfieldSST  could  be  obtained  using  the  full  satellite 
swath  with  a  30%  improvement  in  accuracy  at  large  satellite  zenith  angles  and  that  a  fairly  aggressive  pre-filtering 
scheme  could  help  attain  the  desired  accuracy  of  0.3°K  or  better  using  over  75%  of  the  buoy  matchup  data. 

Keywords:  sea  surface  temperature,  VIIRS  satellite,  remote  sensing,  full  swath. 


1.  INTRODUCTION 

The  ultimate  goal  of  the  prediction  of  Sea  Surface  Temperature  (SST)  from  satellite  data  is  to  attain  an  accuracy  of 
0.3°K  or  better  when  compared  to  floating  or  drifting  buoys  located  around  the  globe.  Current  daytime  SST  algorithms 
are  able  to  routinely  achieve  an  accuracy  of  0.5°K  for  satellite  zenith  angles  up  to  53°.  The  full  scan  swath  of  VIIRS 
(Visible  Infrared  Imaging  Radiometer  Suite)  results  in  satellite  zenith  angles  up  to  70°,  so  that  successful  retrieval  of 
SST  from  VIIRS  at  these  higher  satellite  zenith  angles  would  greatly  increase  global  coverage.  However,  the  accuracy 
of  the  best  SST  algorithm  steadily  degrades  to  nearly  0.7°K  as  the  satellite  zenith  angle  reaches  its  upper  limit,  due 
mostly  to  the  effects  of  increased  atmospheric  path  length. 

Ideally,  the  goal  of  0.3 °K  accuracy  for  angles  up  to  70°  is  to  be  achieved  by  using  only  data  retrieved  from  satellites 
such  as  VIIRS.  Historically,  the  MCSST  algorithm  was  derived  with  this  approach  in  mind  [1],  where  the  coefficient  of 
the  AT=  TirT12  term  is  known  in  the  literature  as  Gamma.  Here  Til  and  T12  are  the  brightness  temperatures  at  11 
microns  (10.7)  and  12  microns,  respectively.  However,  the  use  of  either  climatology  or  gap-free  first  guess  temperature 
field,  Tf^,  was  invoked  by  Walton  [1],  to  express  Gamma  as  a  linear  function  of  surface  temperature:  Gamma  = 
Constant  *Tfield.  This  was  justified  by  experimental  data  gathered  by  Walton,  who  shows  that  its  incorporation  into 
MCSST  to  accentuate  the  linear  dependence  of  Gamma  on  surface  temperature,  known  as  NLSST,  results  in  better 
accuracy.  This  slight  improvement,  however,  also  resulted  in  the  mingling  of  both  Celsius  and  Kelvin  temperatures  into 
the  NLSST  algorithm.  Previous  VIIRS  SST  algorithm  studies  showed  that  Tfield  should  be  in  Celsius  rather  than  Kelvin 
for  best  accuracy  in  the  daytime  NLSST,  while  the  reverse  is  true  for  nighttime  NLSST. 

Recently,  Cayula  and  May  [2]  revisited  this  problematic  inconsistency  in  the  choice  of  units  by  investigating  the  effect 
of  adding  an  offset  from  0°  to  300°  to  Tfield  in  both  the  daytime  and  nighttime  NLSST.  This  new  version,  coined  the 
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expanded  NLSST,  shows  a  slightly  better  accuracy  with  offsets  of  16.2°  and  43.3°  for  daytime  and  nighttime  expanded 
NLSST,  respectively. 

The  use  of  Tfieid  by  Walton  [1]  to  characterize  the  Gamma  coefficient  was  in  effect  equivalent  to  using  it  as  a  predictor. 
However,  the  full  information  provided  by  Tfieid  is  somewhat  tamed  by  the  fact  that  it  is  a  multiplier  to  another  predictor, 
the  AT  term,  which  influences  their  synergy  as  a  combined  predictor.  We  were  therefore  motivated  to  investigate  the  use 
of  Tf^id  as  an  additional,  separate  predictor  in  an  effort  to  extract  the  maximum  amount  of  information  it  contains  with  an 
eye  to  improving  the  accuracy  of  SST  predictions  up  to  53°  Sun  zenith  angles,  as  well  as  exploring  the  full  swath  up  to 
70°.  Only  the  MCSST  and  NLSST  daytime  algorithms  were  used  in  this  investigation. 

After  first  writing  out  the  equations  for  MCSST,  NLSST  and  TfieldSST,  the  SST  algorithm  presented  here  that  includes 
Tfieid  as  a  separate  predictor,  our  approach  is  guided  visually  through  gradual  attempts  to  increase  the  agreement  between 
scatterplots  involving  buoy  data,  Tbuoy-Tn  versus  AT=  Tn-Ti2,  and  scatterplots  involving  SST  predictions,  SST-Tn 
versus  AT=  Tn-Ti2.  Here,  Tbuoy  represents  temperatures  in  Kelvin  units  recorded  at  buoy  locations  around  the  globe, 
while  SST  are  predictions  from  a  particular  SST  algorithm  and  are  also  in  Kelvin  units.  Three  characterizations  of  Tfieid 
are  used,  differing  in  their  spatial  resolution,  in  order  to  assess  the  associated  change  in  accuracy  generated  by 
TfieldSST.  Analysis  of  relationships  between  the  regression  coefficients  is  performed,  resulting  in  a  satisfying  and 
intuitive  representation  of  the  TfieldSST  algorithm  in  terms  of  SSTMC  and  Tfieid . 


2.  RESULTS 


We  begin  by  first  presenting  the  equations  for  MCSST,  NLSST  and  TfieldSST  algorithms  to  define  the  notation  and 
conventions  for  the  terminology  used  in  this  paper: 


MC 

4 


ssr«c  =  c,'K:?;1+c2'K(7;1-r12)+c3*Ks(7;1-7'1!)-c. 
sstnl  =  c,“t, , + cf Tr„u  (r„-r12)+ cf-sfr ,  -  t,  2)-  c. 


NL 

4 


SSTT.,u=c^,+c?'u(Tu-Tl2)+cf'uS(Tu-T,2)-clf'u  +cf"dT] 


field 


(1) 

(2) 

(3) 


where  S  =  sec  0zenith  —  1  .  Note  the  negative  sign  in  front  of  the  c4  coefficients,  in  order  to  be  able  to  only  have  positive 
values  in  coefficient  tables  to  be  presented  later. 

A  buoy  data  set  for  the  month  of  June  2012  created  by  NAVOCEANO  [2]  was  used  to  produce  all  of  the  results 
presented  here.  It  consisted  of  115,036  original  daytime  points.  Routine  NAVOCEANO  filtering  standards  resulted  in 
61,782  points  from  0°  to  53°  ,  and  97,496  points  from  0°  to  70°,  representing  53.70%  and  84.75%  of  the  total  number  of 
points.  It  can  be  seen  that  over  31.05%  of  the  buoy  data  points  are  thrown  away  when  rejecting  the  zenith  band  from  53° 
to  70°.  The  resulting  regression  coefficients  for  the  MCSST  algorithm  when  the  buoy  and  SST  prediction  temperatures 
are  in  °Celsius  were  as  follows: 

cfi'  =  1 .009  cf;  =  2.475  cfc=  1.282  cf:  =  274.9  (4) 

while,  when  the  buoy  and  SST  prediction  temperatures  are  in  °Kelvin  (=  °Celsius  +  273.15°): 

cfc  =1.009  cfic  =  2.475  cfc=  1.282  cfc=1.75  (5) 

Note  that  only  the  value  of  c^C  is  affected  when  comparing  MCSST  predictions  against  different  buoy  temperature 
scales.  Subtracting  Tn  from  both  sides  of  this  last  equation,  the  MCSST  algorithm  can  be  rewritten  as: 


Proc.  of  SPIE  Vol.  8724  87240R-2 


Downloaded  From:  http://proceedings.spiedigitallibrary.org/  on  07/11/2013  Terms  of  Use:  http://spiedl.org/terms 


-1>1.  +cfc(7^ , -7’12)+C3"c5'(7'1  ,  -7^2)-c"c 


(6) 


or,  in  a  form  suitable  for  scatterplot  analysis: 

SSTmc -Tn  =  (c2mc  +cfcs\Tn-Tl2)+Yoffset  (7) 


where  Yoffset  =  (c™c  -1  )rn  -c^c  is  the  narrow  range  of  vertical  intercepts  when  TlrTi2  is  equal  to  0.0,  since  it  depends 
on  Tn.  It  is  narrow  because  its  Tn  coefficient  in  our  case  is  0.009.  As  will  be  seen,  the  value  of  this  offset  is  not  as 
important  as  are  the  potential  clues  that  can  be  gathered  from  comparing  the  shape  and  structure  of  the  buoy  data  and 
SST  predictions  scatterplots.  Side-by-side  comparison  should  reveal  the  similarities  and  differences  between  the  buoy 
data  and  SST  predictions,  as  well  as  suggest  an  approach  to  improve  the  accuracy  of  SST  algorithms.  Figure  1  clearly 
displays  the  challenge  of  SST  predictions  of  sea  surface  temperature.  The  filtering  parameters  used  are  shown  in  the 
upper  left  and  the  buoy  data  scatterplot,  Tbuoy-Tn  versus  Tn-Ti2  as  Figure  la.  The  lower  plots,  Figure  lb-c,  show  the 
scatterplots  predicted  by  MCSST  and  NLSST.  The  color  bar  represents  the  absolute  value  of  the  zenith  angle  in  10° 
bands,  and  the  color  plotting  sequence  is  from  blue  to  red.  Data  is  only  shown  for  zenith  angles  from  0°  to  53°  . 


BUOY 


Filtering  Parameters 

IlfieicfBuoy |  <2.0 
|SST-Buoy  |  <4.0 
|  Zenith  |  <  53  deg 
Adistance  <  25  km 
Atime  <  4  hours 


MCSST 


Filtered  Buoy  Data  (Zenith  Color  Code) 


Data  displays 
scatter 

Goal  is  to  capture  the 
features  and  structure  of 

the  distribution. 


NLSST 


Figure  1.  Visual  comparison  between  scatterplots  involving  buoy  data  and  SST  predictions  for  Sun  zenith  angles  up  to  53°. 


The  most  striking  visual  differences  occur  when  comparing  the  buoy  data  scatterplot  with  the  MCSST  scatterplot  on  the 
bottom  left.  The  “pinch  effect”  and  the  divergence  of  the  separate  10°  zenith  bands  can  be  identified  with  the  form  of  the 
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MCSST  equation.  The  main  discrepancy  appears  to  be  the  absence  of  “scatter”  or  randomness  when  compared  to  the 
buoy  data  scatterplot.  Increasing  this  randomness  in  the  scatterplots  of  SST  predictions  will  result  in  lower  rms  error 
when  compared  to  buoy  data.  Interestingly,  the  inclusion  of  Tfieid  in  the  NLSST  equation  as  part  of  the  second  term 
somehow  introduced  a  slight  amount  of  randomness  as  can  be  seen  in  the  lower  right  NLSST  scatterplot.  The  “pinch 
effect”  is  still  there,  but  moved  up  slightly  due  to  the  non-linearity  introduced  by  using  Tfield  as  a  multiplier. 

As  strongly  suggested  by  this  visual  comparison  between  buoy  and  SST  predictions,  efforts  to  reduce  the  rms  errors 
down  to  0.3°K  from  the  present  daytime  values  hovering  around  0.5°K  will  undoubtedly  require  the  introduction  of 
scatter  and  randomness  into  the  scatterplots  of  future  SST  algorithms.  The  question  therefore  arises  as  to  whether,  if  used 
as  a  separate  predictor  in  the  MCSST  algorithm,  Tfieid  can  provide  more  pronounced  scatter  and  randomness  to  reduce 
the  rms  error  between  SST  retrievals  and  buoy  data?  However,  inclusion  of  an  additional  predictor  requires  precaution 
vis-a-vis  issues  of  multi-collinearity,  as  well  as  measures  of  its  accuracy. 

Characterization  of  Tfidd 

For  this  study,  three  characterizations  of  Tfieid  are  used  in  order  to  gauge  the  effects  of  increasing  spatial  resolution  and 
accuracy:  Clim,  K100  and  K10.  Both  Clim  and  K100  have  100  km  spatial  resolution,  while  K10  enjoys  a  10  km 
resolution.  The  global  rms  errors  for  each  characterization  are  taken  here  to  be  roughly  1.0°K,  0.7°K  and  0.55°K, 
respectively  [2].  Both  K100  and  K10  are  updated  daily  at  NAVOCEANO,  where  there  is  an  effort  to  create  K2  at  2  km 
spatial  resolution.  It  should  be  pointed  out  here  that  neither  K100  nor  K10  ingest  any  buoy-derived  SST  data,  so  they  are 
strictly  based  on  previous  satellite-derived  SST  retrievals.  The  K 100  is  updated  at  least  once  a  day  with  SST  retrievals 
produced  since  the  last  update.  The  K10  is  similar  to  the  K100,  but  the  satellite -derived  SST  data  are  weighted  based  on 
estimated  reliability.  All  three  characterizations  should  provide  additional  information  to  the  MCSST  algorithm  to 
reduce  its  present  rms  error  by  adding  Tfieid  as  a  separate  predictor.  As  previously  mentioned,  this  new  version  of 
MCSST  is  referred  to  here  as  TfieldSST.  However,  if  the  additional  information  provided  by  Tfieidis  redundant,  issues  of 
multi-collinearity  will  surface  and  perhaps  make  the  algorithm  unstable. 

Multicollinearity  of  Tfidd  with  Existing  Predictors 

The  fact  that  the  inclusion  of  Tfieid  as  a  product  term  in  NLSST  generated  randomness  in  the  corresponding  scatterplot 
invites  an  investigation  into  considering  its  effect  as  a  separate  predictor.  Before  we  do  this,  however,  issues  of  its 
potential  multicollinearity  must  be  addressed.  Collinearity  is  defined  as  the  correlation  among  the  predictors  in  a 
multiple  regression,  and  it  therefore  entangles  the  effects  of  the  predictors,  complicating  the  interpretation. 

“Multicollinearity  is  a  statistical  phenomenon  in  which  two  or  more  predictor  variables  in  a  multiple  regression  model 
are  highly  correlated,  meaning  that  one  can  be  linearly  predicted  from  the  others  with  a  non-trivial  degree  of  accuracy.  In 
this  situation  the  coefficient  estimates  may  change  erratically  in  response  to  small  changes  in  the  model  or  the  data. 
Multicollinearity  does  not  reduce  the  predictive  power  or  reliability  of  the  model  as  a  whole,  at  least  within  the  sample 
data  themselves;  it  only  affects  calculations  regarding  individual  predictors.  That  is,  a  multiple  regression  model  with 
correlated  predictors  can  indicate  how  well  the  entire  bundle  of  predictors  predicts  the  outcome  variable,  but  it  may  not 
give  valid  results  about  any  individual  predictor,  or  about  which  predictors  are  redundant  with  respect  to  others.  A  high 
degree  of  multicollinearity  can  also  prevent  computer  software  packages  from  performing  the  matrix  inversion  required 
for  computing  the  regression  coefficients,  or  it  may  make  the  results  of  that  inversion  inaccurate”  [3]. 

If  the  correlation  is  too  high,  then  there  is  no  unique  solution  to  the  regression  coefficients.  There  is  no  rigorous 
threshold  as  to  what  “too  high”  means.  However,  unstable  results  are  expected  if  no  additional  information  is  associated 
with  the  new  predictor,  and  existing  coefficients  change  markedly  when  collinear  predictor  is  added.  Nevertheless, 
stability  may  still  occur  if  the  correlation  with  the  dependent  variable  is  greater  than  the  correlation  with  existing 
predictors. 

Figure  2  displays  scatterplots  of  each  of  the  three  characterizations  of  Tfieid  versus  T1L  The  plots  on  the  left  side  are  color- 
coded  with  respect  to  the  absolute  value  of  the  Sun  zenith  angle,  while  the  ones  on  the  right  are  color-coded  with  respect 
to  the  absolute  value  of  latitude.  The  corresponding  correlation  coefficients  are  shown  at  the  far  left.  It  can  be  seen  that  a 
high  degree  of  correlation  exists  for  all  three  cases,  the  highest  being  associated  with  the  greatest  spatial  resolution,  K10. 
Therefore,  inconsistent  predictions  may  be  generated,  if  the  additional  information  provided  by  Tfieid  is  highly  redundant 
compared  to  the  satellite  data. 
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Tfieid  =  Clim 
R  =  0.96603 


Tfieid  =  K100 
R  =  0.96830 


Tfieid  =K10 

R  =  0.96865 


|Zenith|  Color  Code 


|Latitude|  Color  Code 


|Latitude|  Color  Code 


|Latitude|  Color  Code 


Figure  2.  Scatterplots  of  the  three  Tfieid  characterizations  versus  Tn,  and  associated  correlation  coefficients. 


Figure  2  reveals  that  all  the  scatterplots  display  very  high  correlation  between  Tfield  and  Tn.  As  expected  due  to  the 
increased  spatial  resolution  and  accuracy,  the  correlation  coefficient  R  is  very  high  and  slowly  increases  from  Tfield  = 
Clim  to  K100  to  K10.  These  high  values  of  R  should  result  in  competition  between  the  values  of  the  coefficients  of  Tn 

and  Tfield:  cffield  and  values  should  move  in  opposite  directions  in  a  complementary  fashion.  The  larger  issue  is 

whether  or  not  the  additional  information  provided  by  Tfield  is  redundant,  resulting  in  an  unstable  algorithm.  Only  an 
implementation  of  the  algorithm  is  able  to  address  this  issue. 


Effect  of  Gradual  Increase  ofTfidd  Contribution  to  NLSST 

Armed  with  the  above  warning  signs,  we  performed  an  investigation  similar  to  Cayula  and  May  [2]  using  NLSST,  but 
adding  an  offset  Tdiffoffset  to  the  AT  and  SAT  terms  instead  of  Tfield. 

SSTNL=CTu+C2%jTu -T]2+Tdiffoffset)  +  c?LS(Tu -Tl2+Tdiffoffset)-cNAL  (8) 
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Since  Tfield  is  a  multiplier  of  AT,  this  will  allow  us  to  visually  evaluate  the  effect  of  gradually  increasing  the  effect  of 
T fjeid  and  to  gauge  its  effect  on  the  pattern  of  a  scatterplot  of  SSTNL-Tn  versus  TirT12. 


Tdiffoffset  = -1  rmserror  =  0.62069 


Tdiffoffset  =  +5  rmserror  =  0.45413 


Figure  3.  Effects  on  NLSST-generated  scatterplots  due  to  addition  of  Tdiffoffset  to  Tn-T12. 

The  original  NLSST  with  no  Tdiffoffset  in  Figure  3b.  The  upper  left  plot,  Figure  3  a,  shows  the  effect  of  Tdiffoffset  =  -1, 
which  mostly  serves  to  translate  the  “pinch  effect”  horizontally  to  where  TirT12  =  1.  However,  increasing  Tdiffoffset  to 
5  and  10  Figures  3c  and  3d  results  in  increased  scatter  and  randomness,  as  well  as  decreasing  the  rms  error  from  0.52546 
to  0.45413  to  0.44501.  These  are  the  effects  mentioned  previously  that  we  consider  to  be  essential  ingredients  for  future 
efforts  to  reduce  rms  error.  These  desired  trends  are  a  strong  indication  that  Tfield  could  be  the  vehicle  to  introduce  scatter 
and  randomness  towards  this  goal  and  thereby  reduce  the  rms  error  of  future  SST  algorithms. 


Adding  Tfidd  as  a  Separate  Predictor  in  MCSST 

From  this  point  forward,  MCSST  will  be  used  instead  of  NLSST  in  order  to  clearly  isolate  the  effects  of  adding  Tfield  as  a 
separate  predictor.  In  this  section,  the  effects  of  the  level  of  aggressiveness  of  pre-filtering  of  the  buoy  dataset  for  all 
three  characterizations  of  Tfield  will  be  tabulated  and  overall  relationships  between  regression  coefficients  that  emerge 
from  this  investigation  will  be  presented.  In  all  of  the  following  tables,  c4  values  are  shown  for  SST  predictions  against 
buoy  temperatures  in  °Celsius  in  order  to  clearly  see  the  trend  in  its  value  due  to  its  higher  magnitude. 

Table  1  is  a  repository  for  18  runs  of  TfieldSST,  which  includes  Tfieid  as  a  separate  predictor.  These  runs  are  broken 
down  as  follows:  for  each  the  three  characterizations  of  Tfield  (Clim,  K100,  K10),  6  runs  representing  aggressive  to  non- 
aggressive  pre-filtering  of  the  buoy  data  with  respect  to  Tfieid  were  performed.  The  filtering  bands  are  shown  in  the 
leftmost  column  of  Table  1.  The  TfieldSST  equation  is  shown  at  the  top  without  superscripts  on  the  coefficients  due  to 
lack  of  space  in  trying  to  line  them  up  with  their  associated  column.  Each  column  represents  the  values  of  the  TfieldSST 
coefficients.  The  sum  of  cx  and  c5,  which  we  anticipate  to  be  complementary  due  to  the  high  correlation  between  Tn  and 
Tf^id,  is  shown  in  the  next-to-last  column,  while  the  associated  SSTTfield  rms  error  is  displayed  as  the  end  column. 

The  overall  gradual  trends  in  the  values  of  the  coefficients  serve  to  ease  the  initial  concerns  about  unstability  of  the 
TfieldSST  algorithm.  The  rms  error  gradually  increases  as  the  aggressiveness  of  the  pre-filtering  diminishes  for  all  three 


Proc.  of  SPIE  Vol.  8724  87240R-6 


Downloaded  From:  http://proceedings.spiedigitallibrary.org/  on  07/11/2013  Terms  of  Use:  http://spiedl.org/terms 


characterizations  of  Tfield.  It  can  be  seen  that  the  rms  error  also  diminishes  as  Tfieid  increases  in  accuracy,  from  Clim  to 
K100  to  K10.  Remarkably,  in  the  column  highlighted  in  blue  and  a  strong  confirmation  of  the  stability  of  the  TfieldSST 
algorithm,  the  sum  of  Ci  and  c5  appears  to  be  constant,  of  magnitude  very  close  to  1.0,  illustrating  their  complementarity 
due  to  their  high  correlation.  In  what  follows,  all  rms  errors  will  be  understood  to  be  reported  in  Kelvin  units.  The  rows 
corresponding  to  the  value  of  2.0  for  pre-filtering,  presently  being  used,  are  also  highlighted  in  blue. 


One  should  keep  in  mind  that,  although  aggressive  pre-filtering  serves  to  reduce  the  TfieldSST  rms  error,  it  can  also 
severely  reduce  the  number  of  buoy  data  points  with  which  the  regression  analysis  can  be  performed,  resulting  in  poor 
statistics.  A  trade-off  analysis  between  these  competing  factors  is  discussed  later,  where  it  will  be  pointed  out  that  the 
MCSST  algorithm  is  much  less  sensitive  to  aggressive  pre-filtering  in  reducing  its  “estimated”  rms  error.  Pre-filtering 
does  not  change  the  SST  retrievals  that  are  produced  by  MCSST  or  any  other  SST  algorithm;  rather,  it  only  changes  the 
buoy  data  set  on  which  statistics  are  being  estimated.  In  contrast,  it  will  be  shown  that  the  inclusion  of  Tfield  as  a  separate 
predictor  produces  “genuine”  reduction  of  the  TfieldSST  rms  error  with  increased  pre-filtering  aggressiveness. 


Clim 


Filtering 

+  C2  (T\  1 

1 

^I2)+C3^(^11 

1 

^12)  C4  “ 

I 

<  K100 

Clim 

SST  =  cfiu 

field) 

K10 

1 

1 

1 

1 

Ci  +  C5 

rmserror 

|Tfield'Buoy|  <0.1 

0.021 

0.057 

0.026 

5.921 

0.978 

0.999 

0.07175 

|Tfieid-Buoy|  <0.2 

0.074 

0.194 

0.088 

20.33 

0.927 

1.001 

0.12447 

|Tfield-Buoy|  <0.5 

0.245 

0.643 

0.261 

66.82 

0.759 

1.004 

0.25600 

IT^m-BuovI  <1.0 

0.471 

1.218 

0.535 

128.5 

0.541 

1.012 

0.36675 

[  |Tfl^-Buovl  <2.0 

0.649 

1.663 

0.736 

177.0 

0.364 

1.013 

0.45077  I 

|Tfieid-Buoy|  <5.0 

0.752 

1.914 

0.853 

204.9 

0.259 

1.011 

0.49914 

K100 

1 

i 

I 

I 

1 

1 

C1  +  C5 

rmserror 

|Tfjeld-Buoy  |  <0.1 

0.009 

0.018 

0.005 

2.458 

0.991 

1.000 

0.06968 

|Tfieid-Buoy|  <0.2 

0.043 

0.096 

0.055 

11.86 

0.957 

1.000 

0.12447 

1  Tfieid-Buoy  |  <0.5 

0.195 

0.475 

0.229 

53.23 

0.806 

1.001 

0.25554 

iTfi-w-BuovI  <1.0 

0.389 

0.958 

0.467 

106.1 

0.615 

1.004 

0.36104 

1  |TfioW-Buoy|  <2.0 

0.527 

1.310 

0.633 

143.6 

0.479 

1.006 

0.42815  1 

|Tfield-Buoy|  <5.0 

0.585 

1.462 

0.697 

159.5 

0.420 

1.005 

0.46373 

K10 

I 

I 

i 

1 

1 

Ci  +  C5 

rmserror 

|Tfield-Buoy|  <0.1 

0.013 

0.032 

0.011 

3.600 

0.986 

0.999 

0.06991 

1  Fieid-Buoy  |  <0.2 

0.055 

0.130 

0.063 

15.10 

0.945 

1.000 

0.12362 

1  Feid-Buoy  |  <0.5 

0.188 

0.454 

0.239 

51.21 

0.814 

1.002 

0.24467 

|TfiQlH-Buovl  <  1.0 

0.326 

0.810 

0.406 

89.01 

0.677 

1.003 

0.33348 

I  |TfiDlH-Buoy|  <2.0 

0.419 

1.049 

0.518 

114.3 

0.585 

1.004 

0.39268  1 

lFjeld"Buoy  |  <5.0 

0.467 

1.169 

0.571 

127.2 

0.539 

1.006 

0.43082 

Table  1.  Effects  of  pre-filtering  buoy  data  and  increased  spatial  resolution  of  Tfield 
using  MCSST  for  Sun  zenith  angles  from  0°  to  53°. 

Changes  in  Scatterplots  due  to  Tfidd  as  a  Separate  Predictor 

The  addition  of  Tfieid  as  a  separate  predictor  also  results  in  the  addition  of  significant  randomness  to  the  scatterplots 
previously  shown  in  Figure  1,  a  crucial  feature  mentioned  earlier  towards  increasing  the  agreement  between  SST 
predictions  and  the  buoy  data.  Figure  4  compares  the  buoy  data  scatterplot  versus  the  scatterplot  predicted  by  MCSST, 
basically  TfieldSST  without  Tfield  as  a  separate  predictor,  for  zenith  angles  from  0°  to  70°,  in  10°  bands. 

Because  the  color  plotting  routine  of  each  band  occurs  sequentially,  some  of  the  scatter  in  one  band  overwrites  the 
previous  band.  As  a  way  for  the  reader  to  assess  the  amount  of  overlap  between  zenith  bands,  the  plotting  sequence  is 
shown  for  both  Blue  to  Red  and  Red  to  Blue  (Figure  4a-b  for  buoy  data,  Figure  4c-d  for  MCSST).  The  resulting  dual 
scatterplots  for  MCSST  reveal  that  there  is  practically  no  overlap  of  zenith  band  scatter,  except  for  low  values  of  TirTi2. 
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In  sharp  contrast,  however,  the  buoy  data  scatterplot  shows  pronounced  overlap  of  zenith  band  scatter.  Figure  5  is  similar 
to  Figure  4,  but  with  the  lower  scatterplots  predicted  by  TfieldSST  using  Tfieid  =  K10,  instead  of  MCSST.  As  can  be 
seen,  features  very  similar  to  the  buoy  data  scatterplot  emerge,  thereby  establishing  a  strong  link  between  reduced  rms 
error  and  the  presence  of  scatter  and  randomness  in  the  TfieldSST-predicted  scatterplot,  as  was  conjectured  earlier. 


Blue-*  Red  COLOR  PLOTTING  Red  -*  Blue  COLOR  PLOTTING 


Figure  4.  Comparison  of  MCSST  predictions  with  Buoy  Data  scatterplots,  without  Tfield, 
for  Sun  zenith  angles  from  0°  to  70°. 


Blue-*  Red  COLOR  PLOTTING  Red  -*  Blue  COLOR  PLOTTING 


Figure  5.  Comparison  of  MCSST  predictions  with  Buoy  Data  scatterplots,  with  Tfield  =  K10, 
for  Sun  zenith  angles  from  0°  to  70°. 
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Dependence  of  Regression  Coefficients  on  Zenith  and  Latitude 

The  extreme  stability  of  the  TfieldSST  algorithm  when  Tfield  =  K10  is  also  evident  from  the  gradual  changes  in  the  values 
of  the  regression  coefficients  as  either  zenith  angle  or  latitude  is  broken  up  into  equally  spaced  bands.  Besides  the  steady 
increase  of  rms  error  as  the  zenith  angle  increases,  Table  2  below  shows  this  general  gradual  change  in  the  values  of  the 
coefficients  with  a  single  glaring  exception,  such  as  the  value  of  -0.011  for  c3  for  the  20°-30°  zenith  band.  It  is  possible 
that  some  outlier  buoy  data  values  creeped  into  the  original  creation  of  the  data  file. 

In  addition,  global  coefficient  values  over  the  entire  zenith  span  from  0°  to  70°  are  shown  in  order  to  compare  them  with 
the  global  coefficient  values  derived  from  MCSST  shown  in  green.  Again,  the  rms  error  shows  a  gradual  increase  in  rms 
error  over  the  10°  zenith  bands  when  using  these  global  coefficients,  contrasted  with  the  MCSST  rms  errors  in  green. 


From  GUI  results 


|  Zenith  |  Band 

#pts 

rmserror 

bias 

Cl 

c2  c3 

c4  c5 

(Ci+C5) 

+60  to  +70 

21213 

0.41834 

-8.37e-13 

0.230 

0.292  0.296 

62.7  0.794 

1.024 

+50  to  +60 

19667 

0.39780 

-2.03e-13 

0.296 

0.621  0.349 

80.9  0.721 

1.017 

+40  to  +50 

15353 

0.39463 

2.83e-13 

0.343 

0.814  0.447 

93.4  0.665 

1.008 

+30  to  +40 

12415 

0.38942 

-9.51e-13 

0.407 

1.004  0.742 

110.8  0.594 

1.001 

+20  to  +30 

10231 

0.38372 

-2.92e-12 

0.469 

1.264  -0.011 

127.7  0.533 

1.002 

+10  to  +20 

9477 

0.37942 

7.47e-13 

0.521 

1.325  1.087 

141.9  0.482 

1.003 

0  to +10 

9151 

0.37119 

-6.32e-13 

0.493 

1.292  0.922 

134.3  0.506 

0.999 
^  _ 

Using  global  coefficients 

97507 

0.41055 

3.66e-12 

0.294 

0.610  0.388 

80.10  0.719 

1.013 

no  K10 

0.65727 

1.09e-ll 

1.034 

2.204  1.363 

281.8  0.000 

1.034 

|  Zenith  |  Band 

#pts 

rmserror 

bias 

rmserror 

bias 

+60  to  +70 

21213 

0.44234 

-0.04142 

0.86911 

-0.16031 

+50  to  +60 

19667 

0.39847 

-0.00778 

0.63142 

0.01057 

+40  to  +50 

15353 

0.39828 

0.00676 

0.58992 

0.05601 

+30  to  +40 

12415 

0.40274 

0.02223 

0.58093 

0.06392 

+20  to  +30 

10231 

0.40334 

0.01332 

0.53541 

0.05429 

+10  to  +20 

9477 

0.40560 

0.03403 

0.50201 

0.06380 

0  to  +10 

9151 

0.39493 

0.02127 

0.52607 

0.04228 

Table  2.  Dependence  of  regression  coefficients  on  Sun  zenith  angle,  with  Tfield  =  K10. 

The  same  kind  of  analysis  was  performed  in  20°  latitude  bands  and  the  values  of  the  resulting  coefficients  are  presented 
in  Table  3.  A  gradual  increase  in  rms  error  can  be  seen  as  the  latitude  bands  are  further  removed  from  the  equator,  except 
for  the  two  extreme  bands  near  the  poles.  The  sum  of  Ci  and  c5  again  shows  remarkable  constancy,  except  for  the  value 
of  0.914  at  the  equatorial  band. 

As  was  done  in  Table  2,  the  values  of  global  coefficients  for  both  TfieldSST  with  K10  and  MCSST  (green)  are  displayed 
for  comparison.  The  use  of  these  global  coefficients  to  calculate  the  rms  errors  over  the  latitude  bands  reveals  several 
outliers,  whose  large  biases  are  highlighted  in  blue.  The  rms  errors  associated  with  these  large  biases,  although  still 
within  the  range  of  values  of  the  other  rms  errors,  are  also  relatively  high.  The  MCSST  global  coefficients  also  predicted 
large  biases  (green)  highlighted  in  blue. 
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From  GUI  results 


- \ 


Latitude  Band 

#pts 

rmserror 

bias 

Cl 

c2 

c3 

c4 

C5 

(Cj+Cg) 

+50  to  +70 

7112 

0.36861 

-1.38e-12 

0.375 

0.980 

0.314 

102.2 

0.626 

1.001 

+30  to  +50 

17404 

0.39976 

5.79e-13 

0.513 

1.191 

0.626 

139.8 

0.497 

1.010 

+10  to  +30 

12910 

0.3532 

-6.44e-13 

0.279 

0.672 

0.339 

76.3 

0.731 

1.010 

-10  to  +10 

4856 

0.3521 

-1.66e-13 

0.242 

0.795 

0.240 

63.9 

0.672 

0.914 

-30  to  -10 

8947 

0.34733 

-1.81e-13 

0.432 

1.080 

0.554 

117.9 

0.580 

1.012 

-50  to  -30 

9098 

0.38322 

-6.67e-13 

0.514 

1.054 

0.612 

140.0 

0.501 

1.015 

-70  to  -50 

1468 

0.36599 

-2.07e-14 

0.359 

0.777 

0.359 

98.1 

0.673 

1.032 
k  J 

Using  global  coefficients 

61795 

0.3894 

3.05e-12 

0.416 

1.039 

0.515 

113.4 

0.589 

1.005 

no  K10 

0.5305 

6.62e-12 

1.006 

2.550 

1.174 

274.1 

0.000 

1.006 

Latitude  Band 

#pts 

rmserror 

bias 

rmserror 

bias 

+50  to  +70 

7112 

0.37013 

-0.02832 

0.5056 

-0.07924 

+30  to  +50 

17404 

0.40451 

10.12566  I 

0.51462 

'0.11358" 

+10  to  +30 

12910 

0.36430 

-0.01473 

0.57448 

-0.11061 

-10  to  +10 

4856 

0.38496 

Umsss  J 

0.58587 

,-0.11294, 

-30  to  -10 

8947 

0.34863 

-0.07716 

0.47300 

0.00595 

-50  to  -30 

9098 

0.38908 

-0.06359 

0.49011 

0.05780 

-70  to  -50 

1468 

0.37408 

0.03549 

0.47590 

-0.00664 

Table  3.  Dependence  of  regression  coefficients  on  latitude  bands  for  Sun  zenith  angles  from  0°  to  53°,  with  Tfield  =  K10. 


Reduction  ofrms  error  as  a  function  of  Sun  zenith  angle  bands 

Present  SST  algorithms  only  consider  Sun  zenith  angles  from  0°  to  53°.  In  an  effort  to  quantify  the  improvement  in 
accuracy  afforded  with  the  use  of  Tfield  as  a  separate  predictor,  Table  4  displays  the  rms  errors  predicted  by  the  MCSST 
and  the  TfieldSST  algorithms  when  the  three  characterizations  of  Tfield  are  used.  The  percent  reduction  for  each 
characterization  is  shown  in  green.  The  reduction  in  rms  error  in  the  zenith  band  from  0°  to  53°  is  clearly  seen  from  left 
to  right,  as  the  spatial  resolution  of  Tfield  increases  from  Clim  to  K100  to  K10.  In  addition,  results  for  the  zenith  bands  of 
53°  to  70°  and  0°  to  70°  representing  the  missing  part  of  the  full  swath  are  shown  in  the  last  2  rows  of  Table  4.  The  rms 
error  percent  reduction  over  the  53°  to  70°  zenith  band  is  quite  remarkable  as  the  rms  error  drops  from  0.73712  down  to 
0.41298  for  Tfield  =  K10.  The  full  swath  zenith  band  of  0°  to  70°  shows  a  percent  reduction  of  over  30%  with  K100  and 
even  more  with  K10. 


Zenith 

MCSST 

Tfield  =  Clim 

Tfield  = 

K100 

Tfield  = 

K10 

0°  to  53° 

0.53052 

0.45077  15.0% 

0.42551 

19.8% 

0.38937 

26.6% 

53°  to  70° 

0.73712 

0.51667  29.9% 

0.45622 

38.1% 

0.41298 

43.9% 

0°  to  70° 

0.65627 

0.49989  23.8% 

0.45507 

30.7% 

0.41053 

37.4% 

Table  4.  Percent  reduction  in  rms  error  (green)  over  Sun  zenith  angle  bands  for  the  three  characterizations  of  Tfield. 
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Global  Improvement  of  SST predictions  against  Buoy  data 

Another  way  to  visualize  these  genuine  improvements  over  the  3  zenith  bands  in  Table  4  is  to  produce  worldmaps  of 
SSTTfieid-Tbuoy  corresponding  to  these  zenith  bands  side-by-side  with  worldmaps  of  SSTMC-Tbuoy.  Figures  6,  7  and  8  use 
the  Blue  to  Red  and  Red  to  Blue  method  to  plot  results  for  0°  to  53°,  53°  to  70°  and  0°  to  70°,  respectively. 


Blue  Red  Color  Plotting 

BUOY  LOCATIONS  (SSTBuoy  Color  Code) 


Red  4  Blue  Color  Plotting 
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Figure  6.  Global  map  of  SSTMC-Tbuoy  (top)  and  SSTTfield-Tbuoy  (bottom)  for  Sun  zenith  angles  from  0°  to  53°. 
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Figure  7.  Global  map  of  SSTMC-Tbuoy  (top)  and  SSTxfield-Tbuoy  (bottom)  for  Sun  zenith  angles  from  53°  to  70°. 
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For  all  three  figures,  the  top  two  worldmaps  represent  predictions  from  MCSST,  while  the  bottom  two  worldmaps  are 
the  TfieldSST  predictions  using  Tfieid  =  K10.  The  genuine  reduction  in  rms  errors  described  previously  is  evidenced  by 
the  ability  of  the  TfielSST  algorithm  to  remove  extreme  differences  that  are  present  in  MCSST  results,  both  positive 
(towards  the  red)  and  negative  (towards  the  blue).  Note  the  impressive  confinement  of  differences  to  the  yellow  and  pale 
blue  colors  upon  close  inspection  of  the  TfieldSST  worldmaps. 
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Figure  8.  Global  map  of  SSTMC-Tbuoy  (top)  and  SSTTfield-Tbuoy  (bottom)  for  Sun  zenith  angles  from  0°  to  70°. 


3.  DISCUSSION 


The  previous  results  from  using  Tfieid  as  a  separate  predictor  revealed  a  few  salient  features  and  relationships  that  bear 
discussing.  We  offer  an  interpretation  of  the  results  that  clearly  shows  how  and  why  the  synergy  in  the  combination  of 
daily  satellite  data  with  Tfieid,  a  no-gap  temperature  field  derived  from  previous  satellite  data,  is  capable  of  genuine 
reduction  in  the  rms  error  of  SST  predictions  of  sea  surface  temperature.  The  major  differences  between  pre -filtering  of 
the  buoy  data  set  with  respect  to  Tfield  and  the  use  of  higher  spatial  resolution  versions  of  Tfield  are  highlighted. 
Advantages  of  using  Tfield  as  a  separate  predictor  are  discussed. 


Interpretation  of  Results 

Equations  (1)  and  (3)  for  the  two  SST  algorithms  compared  in  this  paper,  MCSST  and  TfieldSST,  are  again  displayed 
here  for  easy  reference  in  this  discussion: 


SSTUC  =c"cTn+c"c(Tu-Tl2)+c"cS(Tn-T,1)-c. 


MC 

4 


SSTnu=c^,+cY'eU(Tn-Tn)+cf’uS(Tu-TI2)-cY‘-u+cY^T. 


field  ' 


(1) 

(3) 


As  shown  in  Results  section,  the  sum  of  the  Tn  and  Tfieid  coefficients  remained  remarkably  steady  around  the  value  of 
unity: 
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(9) 


Tfield  .  Tfield 

C\  +  c/ 


=  1. 


In  a  search  for  other  relationships  that  could  possibly  offer  additional  insight  into  an  interpretation  of  results,  a  separate 
investigation  of  possible  relationships  between  the  ratios  of  TfieldSST  coefficients  and  MCSST  coefficients  was 
performed.  The  results  are  shown  in  Tables  5  and  6,  which  differ  only  in  the  range  of  the  Sun  zenith  angles.  The  values 
in  red  were  obtained  by  dividing  the  TfieldSST  coefficients  2,  3  and  4  by  the  first  TfieldSST  coefficient,  cffield  . 


cl 

c2 

c3 

c4 

c5 

rmserror 

MCSST 

1.006 

2.541 

1.173 

274.10 

— 

0.5058 

TfieldSST 

0.251 

0.617 

0.312 

68.42 

0.752 

0.2869 

1.000 

2.467 

1.243 

272.58 

— 

Table  5.  Comparison  of  coefficients  of  MCSST  and  TfieldSST,  as  well  as  their  ratios  . 
Sun  zenith  angles  are  from  0°  to  53°. 


cl  c2  c3  c4  c5  rmserror 
MCSST  1.034  2.201  1.360  281.98  —  0.6382 

TfieldSST  0.163  0.326  0.221  44.48  0.844  0.2987 

1.000  2.001  1.356  272.88  — 


Table  6.  Same  as  Table  5,  except  that  Sun  zenith  angles  are  from  0°  to  70° 


Because  of  the  general  trends  in  coefficient  values  presented  in  tabular  form  earlier  and  the  relatively  good  agreement  in 
values  between  the  MCSST  coefficients  and  these  ratios,  the  following  relationship  suggests  itself: 


c. 


MC  _  Tfield  I  Tfield 


/  C| 


or 


fTfield 


=  C 


Tfield  MC 


C, 


Substituting  this  relationship  into  the  equation  for  SSTTfieid: 

CCT  _  Tfield  \MCrp  ,  MC  / y-  _rp  \\nMCc(rr  _rT  nMC}\  nTfield  T 
1  Tfield  ~  C1  Lci  ill"hC2  Vll  I12/~hC3  °V*11  112j  C4  J“hC5  1J 


field 


resulting  in  a  relationship  that  clearly  illustrates  a  linear  weighted  average  between  only  two  predictors: 

_  Tfield  cc^r  .  U  _  Tfield 
007  Tfield  ~L\  007  MC  V1  C1  r  fl 


L  field  • 


(10) 


(11) 


(12) 


The  latter  equation  can  be  rewritten  as  a  linear  function  of  the  difference  between  SSTMC  and  Tfi 


SSTV,U  =Tm  +C^(SST„C -Tf,u).  (13) 

A  very  simple  and  insightful  interpretation  now  emerges.  The  TfieldSST  algorithm  is  searching  for  a  correction  to  the 
background  (a  priori  satellite  data  through  Tfieid)  which  is  a  linear  function  of  the  difference  between  observation  (new 
satellite  data  through  SSTMC)  and  background  (Tfieid).  Of  great  relevance  and  guidance  here  is  the  established  theory 
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about  the  relationship  between  the  individual  variances  for  this  type  of  equation.  If,  for  the  sake  of  simplicity  in  the 
resulting  relationships  for  discussion  purposes,  the  observation  and  background  errors  are  now  assumed  to  be 
uncorrelated  [5]: 


Tfield 

C1 


(14) 


The  optimal  cffield  that  minimizes  the  error  variance  G 2SSJ ^  is  obtained  by  setting  the  derivative  with  respect  to  c^ldd 
equal  to  0.0,  resulting  in: 


J field 


' Tfield 


and 


1  _  rTfield  —  _ 

1  Cj  —  C5  - 


Tfield 


' Tfield 


+  <J\ 


sstm. 


sstM( 


'  Tfield 


+  (J 


sstM( 


(15) 


When  this  optimal  value  for  c^ield  is  substituted  in  Equation  (14),  a  simple  relationship  between  the  three  variances 
emerges: 


sstM( 


2  2 
\G Tfield  +CJSSTmc  J 


G Tfield  + 


7 Tfield 


2  2 
\° Tfield  +<JSSTmc  J 


'SSTm 


or 


1 


1  1 

T2 - +  ~2 - 

G Tfield  GSSTmc 


(16) 


Note  that  the  resulting  variance  of  SSTTfieid  will  always  be  smaller  than  either  variance,  due  to  the  “resistors  in  parallel” 
characteristic  of  this  last  relationship.  Clearly,  should  the  variances  of  Tfieid  and  SSTMC  be  equal,  the  resulting  variance  of 
SSTTfieid  will  be  half  of  each,  resulting  in  its  rms  error  being  a  factor  of  V2 /2  smaller  than  their  rms  errors.  For  example, 
should  the  rms  errors  of  Tfield  and  SSTMC  both  be  equal  to  0.5,  the  resulting  rms  error  of  SSTTfieid  would  be  0.3535,  a 
substantial  reduction  that  is  in  the  neighborhood  of  TfieldSST  rms  errors  reported  in  this  article.  This  interpretation  and 
its  consequences  on  the  variance  of  SSTTfield  correctly  and  neatly  encapsulates  the  nature,  magnitude  and  overall  trend  of 
the  results  presented  earlier. 


A  couple  of  other  interesting  relationships  which  relate  the  rms  errors  between  SSTTfieid  and  its  predictors  SSTMC  and 
Tfieid  immediately  follow  from  the  above  equations: 


_  Tfield  2 


SSTrfa 


=  Ci 


'SSTM( 


SSTn. 


=  ,R 


Tfield 


'SSTM( 


(17) 


and 


2  _  Tfield  2 

aSSTTfield  ~  C5  GTfield 


° Tfield  • 


(18) 


From  this,  a  ballpark  estimate  of  the  percent  reduction  in  SSTMC  rms  error  is  conveniently  given  by  100^1  -  sjclfield  j. 

Introduction  of  a  small  correlation  between  SSTMC  and  Tfieid  errors  with  respect  to  Tbuoy  will  introduce  minimal  changes 
in  these  ballpark  estimates. 
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The  global  rms  error  for  the  Clim,  K100  and  K10  fields  were  estimated  to  be  about  1.0,  0.7  and  0.55,  respectively  [2]. 
Substituting  the  rms  error  for  SSTMC  from  Table  5,  we  arrive  at  SSTTfield  rms  error  values  of  0.450,  0.409  and  0.372  for 
zenith  angles  from  0°  to  53°  .Using  the  SSTMC  rms  error  value  from  Table  6,  SSTTfield  rms  error  values  of  0.538,  0.472 
and  0.4166  are  obtained  when  zenith  angles  from  0°  to  70°  are  considered.  These  values  again  are  consistent  with  the 
results  concerning  the  three  different  characterizations  of  Tfield  were  used:  Clim,  K100  and  K10. 

With  the  above  explanation  of  the  reason  why  the  TfieldSST  rms  error  turns  out  to  be  substantially  smaller  than  either 
the  SSTmc  or  Tfieid  rms  errors,  we  now  turn  to  a  discussion  of  the  achievability  of  the  0.3 °K  threshold  for  rms  error  using 
any  SST  algorithm.  Figure  9  displays  several  curves  which  represent  the  Desired  rms  error  as  a  function  of  both  SSTMC 
and  Tfieid  rms  errors,  using  Equation  (16).  The  bold  red  line  represents  the  Desired  0.3°K  threshold  for  rms  error.  The 
black  circle  corresponds  to  a  Desired  rms  error  of  0.3535  just  discussed  in  the  example  above. 


Figure  9.  Curves  of  Desired  rms  errors. 


The  diagonal  dashed  line  represents  the  path  taken  to  reduce  the  Desired  rms  error  with  a  common  reduction  factor  for 
both  SSTmc  and  Tfield.  Its  intersection  with  the  0.3°K  threshold  curve  (red)  requires  that  each  rms  error  be  reduced  to  a 
value  of  0.3  /  ( V2/2)  or  0.425,  a  reduction  of  about  15%  from  existing  rms  errors  of  roughly  0.5  for  both  SSTMC  and 
Tfieid-  From  the  relationship  between  variances  given  earlier  in  Equation  (16),  the  path  along  the  vertical  dashed  line 
would  require  a  reduction  in  the  SSTMC  rms  error  to  a  value  of  0.375,  while  keeping  the  Tfieid  rms  error  at  0.5.  Due  to  the 
symmetric  nature  of  the  variances  relationship  described  by  Equation  (16)  and  also  illustrated  in  Figure  9,  the  path  along 
the  horizontal  dashed  line  would  require  a  reduction  of  Tfieid  rms  error  to  a  value  of  about  0.375  as  well,  while  keeping 
the  SSTmc  rms  error  at  0.5.  This  represents  a  reduction  of  about  25%  from  existing  Tfieid  rms  error  of  roughly  0.5. 

The  latter  path  seems  to  be  the  easiest  and  quickest  way  to  achieve  the  0.3  °K  threshold,  as  efforts  are  now  under  way  at 
NAVOCEANO  [2]  to  develop  a  K2  characterization  of  Tfield,  with  a  spatial  resolution  of  2  km.  A  decrease  in  Tfield  rms 
error  is  expected  with  this  new  product,  due  to  its  higher  spatial  resolution. 

Interestingly,  if  one  is  willing  to  sacrifice  a  small  percentage  of  the  filtered  buoy  data  points,  the  0.3°K  threshold  can  be 
easily  achieved  presently  without  any  improvements  to  the  present  SSTMc  or  Tfield  rms  errors.  More  aggressive  pre¬ 
filtering  of  the  buoy  data  with  respect  to  Tfield  is  analogous  to  artificially  increasing  the  accuracy  of  Tfield,  at  the  cost  of 
eliminating  some  buoy  data  points.  We  therefore  now  present  the  results  of  a  tradeoff  study  between  pre-filtering 
aggressiveness,  spatial  resolution  of  Tfieid,  range  of  zenith  angles  and  resulting  number  of  filtered  buoy  data  points 
available  for  ingestion  into  the  TfieldSST  algorithm,  with  the  ultimate  goal  of  achieving  the  0.3 °K  threshold  in  Figures 
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10,  11  and  12,  where  runs  were  made  for  maximum  values  of  ITfield  -  Tbuoyl  =  0.1,  0.2,  0.5,  1.0,  2.0,  3.0,  4.0  and  5.0. 
These  figures  only  differ  with  respect  to  zenith  bands:  0°  to  53°,  53°  to  70°  and  full  swath  0°  to  70°. 

For  each  figure,  the  horizontal  axis  represents  the  aggressiveness  of  the  pre -filtering,  increasing  from  right  to  left.  The 
plot  on  the  left  is  the  rms  error  resulting  from  SSTMC  (bold  red),  and  the  three  characterizations  of  Tfieid.  The  plot  on  the 
right  displays  two  variables:  the  percent  of  buoy  data  left  over  after  the  pre-filtering  (solid  lines),  and  the  percent 
contribution  that  Tfield  makes  to  the  TfieldSST  algorithm,  expressed  as  cfield  in  percent  units.  The  0.3°K  threshold  is 

represented  by  the  horizontal  dashed  black  line.  The  vertical  dashed  line  is  a  visual  estimate  of  the  maximum  value  of 
ITfieid  -  Tbuoyl  required  to  achieve  this  goal,  and  is  drawn  for  a  value  of  0.7°K  on  both  plots.  A  noteworthy  feature  is  that 
the  K100  and  K10  characterizations  of  Tfield  seem  to  converge  at  that  value,  while  Clim  requires  more  aggressive 
filtering  to  reach  the  0.3°K  threshold. 


0  <  Zenith  Angles  <  53  0  <  Zenith  Angles  <  53 


Figure  10.  Plots  of  rms  error,  %  of  buoy  data  points  filtered  and  %  of  Tfield  contribution  to  TfieldSST,  (100  c^ldd ), 

for  Sun  zenith  angles  from  0°  to  53°. 


It  will  be  recalled  from  the  tabular  data  presented  in  the  Results  section,  that  overall  trends  showed  that  the  value  of 
cfeld  depended  on  the  aggressiveness  of  the  pre-filtering  of  the  buoy  dataset  against  Tfield,  the  resolution  of  the  Tfieid  data 

(Clim,  K100,  K10),  as  well  as  the  Sun  zenith  angle.  In  general,  for  a  given  level  of  pre-filtering,  the  value  of  cfield 

increased  and  the  rms  error  decreases  as  the  spatial  resolution  of  the  Tfield  data  increased  from  Clim  to  K100  to  K10. 
This  indicates  that  the  TfieldSST  algorithm  automatically  adjusts  the  mix  of  satellite  and  Tfieid  data  to  produce  rather 
drastic,  but  genuine,  increases  in  accuracy  of  about  30%  (equivalent  to  the  multiplicative  factor  of  V2/2  found 
previously).  Results  also  revealed  that  the  value  of  cfldd  increases  with  zenith  angle,  showing  that  present  regression 

algorithms  that  use  only  daily  satellite  data  are,  by  themselves,  less  capable  and  less  accurate  as  the  zenith  angle 
increases.  The  addition  of  a  priori  satellite-derived  data  offered  by  Tfield  generates  meaningful  reductions  in  rms  error. 
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53  <  Zenith  Angles  <  70  53  <  Zenith  Angles  <  70 


Figure  11.  Plots  of  rms  error,  %  of  buoy  data  points  filtered  and  %  of  Tfield  contribution  to  TfieldSST,  (100  cjield ) 

for  Sun  zenith  angles  from  53°  to  70°. 


0  <  Zenith  Angles  <  70  0  <  Zenith  Angles  <  70 


Figure  12.  Plots  of  rms  error,  %  of  buoy  data  points  filtered  and  %  of  Tfield  contribution  to  TfieldSST,  (100  ) 

for  Sun  zenith  angles  from  0°  to  70°. 
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All  of  these  gradual  trends  are  reflected  in  Figures  10  through  12  and  therefore  reinforce  our  assumption  that  TfieldSST 
is  basically  a  linear  weighting  of  SSTMC  (daily  satellite  data)  and  Tfieid  (a  priori  satellite  data).  It  is  only  for  high  values  of 
ITfieid  -  Tbuoyl,  corresponding  to  less  aggressive  pre-filtering,  that  differences  in  SSTTfield  rms  errors  reveal  themselves 
between  each  characterization  of  Tfield.  This  can  be  seen  on  all  the  plots  of  rms  error  vs  ITfield  -  Tbuoyl  in  Figures  10 
through  12.  Unlike  the  Tfieid  rms  errors,  the  SSTMC  rms  error  shown  in  bold  red  is  relatively  insensitive  to  the 
aggressiveness  in  pre-filtering,  for  all  zenith  bands  shown  in  Figures  10  through  12.  For  the  53°  to  70°  band  of  zenith 
angles  shown  in  Figure  11,  the  rms  error  averages  about  0.725,  for  a  combined  full  swath  rms  error  average  of  about 
0.650  shown  in  Figure  12.  These  trends  are  consistent  with  our  earlier  remarks  that  pre-filtering  of  the  buoy  data  with 
respect  to  Tfield  does  not  change  the  SST  retrievals  that  are  produced  by  MCSST  or  any  other  SST  algorithm.  In  sharp 
contrast,  however,  the  gradual  and  real  decrease  in  SSTTfield  rms  error  seen  from  right  to  left  for  all  three 
characterizations  of  Tfield  is  consistent  with  the  artificial  reduction  of  Tfieid  rms  error  achieved  through  gradually  more 
aggressive  pre-filtering. 

The  plots  on  the  right  of  each  figure  are  more  pertinent  to  our  tradeoff  analysis,  however,  where  the  same  vertical  dashed 
line  with  a  maximum  value  of  0.7°K  for  ITfield  -  Tbuoyl  required  to  achieve  the  goal  of  0.3°K  threshold  is  replicated.  It 
should  be  mentioned  that  Clim  is  also  used  for  the  pre-filtering  of  buoy  data  when  using  MCSST,  explaining  why  there 
is  no  solid  red  curve  and  that  the  green  solid  curve  applies  for  MCSST  pre-filtering.  In  addition,  since  there  is  no  Tfield 
term  in  MCSST,  the  red  dashed  curve  for  MCSST  represents  a  0.0  value  for  100  cfield  . 

Figure  10  for  the  0°  to  53°  zenith  band  shows  that  only  about  55%  of  the  data  buoy  points  filtered  through  at  the  less 
aggressive  range  of  ITfield  -  Tbuoyl  for  all  characterizations  of  Tfieid.  In  order  for  the  TfieldSST  algorithm  to  achieve  the 
0.3°K  threshold  required  with  value  0.7°K  of  ITfield  -  Tbuoyl,  only  about  25%  of  the  buoy  data  points  for  Clim  are  left  to 
be  ingested  into  the  algorithm,  while  both  K100  and  K10  only  suffer  a  more  modest  decline  to  42%  and  46%, 
respectively. 

Figure  11  for  the  53°  to  70°  zenith  band  shows  that  only  a  little  over  30%  of  the  data  buoy  points  filtered  through  at  the 
less  aggressive  range  of  ITfield  -  Tbuoyl  for  all  characterizations  of  Tfieid.  In  order  for  the  TfieldSST  algorithm  to  achieve 
the  0.3°K  threshold  required  with  a  maximum  value  of  0.7°K  for  ITfield  -  Tbuoyl,  only  about  15%  of  the  buoy  data  points 
for  Clim  are  left  to  be  ingested  into  the  algorithm,  while  both  K100  and  K10  only  suffer  a  more  modest  decline  to  26% 
and  28%,  respectively. 

Finally,  Figure  12  for  the  0°  to  70°  zenith  band  shows  that  a  little  over  85%  of  the  data  buoy  points  filtered  through  at 
the  less  aggressive  range  of  ITfieid  -  Tbuoyl  for  all  characterizations  of  Tfieid.  In  order  for  the  TfieldSST  algorithm  to 
achieve  the  0.3 °K  threshold  required  with  value  0.7°K  of  ITfieid  -  Tbuoyl  I,  only  about  45%  of  the  buoy  data  points  for 
Clim  are  left  to  be  ingested  into  the  algorithm,  while  both  K100  and  K10  only  suffer  a  more  modest  decline  to  68%  and 
74%,  respectively. 

Summarizing,  for  the  present  0°  to  53°  range  of  zenith  angles  considered  in  SST  algorithms,  a  reduction  of  about  only 
9%  in  filtered  buoy  data  points  results  in  a  drop  in  rms  error  from  0.5  for  MCSST  to  the  goal  of  0.3°K  when  Tfield  =  K10 
is  used  in  the  TfieldSST  algorithm.  Similarly,  for  the  53°  to  70°  range  of  zenith  angles  presently  ignored  in  SST 
algorithms,  a  reduction  of  about  only  2%  in  filtered  buoy  data  points  results  in  a  drop  in  rms  error  from  0.75  for  MCSST 
to  the  goal  of  0.3°K  for  SSTTfield  using  Tfield  =  K10.  Finally,  for  the  full  swath  0°  to  70°  range  of  zenith  angles,  a 
reduction  of  about  only  1 1%  in  filtered  buoy  data  points  results  in  a  drop  in  rms  error  from  0.65  for  SSTMC  to  the  goal  of 
0.3°K  for  SSTTfieid  using  T fieid  =  K10. 

It  should  be  emphasized  at  this  point  that  only  about  55%  of  the  buoy  data  points  are  being  presently  being  used  to 
achieve  a  rms  error  of  about  0.5  for  MCSST  considering  only  the  0°  to  53°  range  of  zenith  angles,  while  use  of 
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TfieldSST  now  provides  about  75%  of  the  buoy  data  points  over  the  full  swath  from  0°  to  70°  range  of  zenith  angles 
with  an  overall  rms  error  of  0.3°K,  attaining  the  seemingly  elusive  goal  of  SST  algorithm  accuracy. 

Our  assumption  about  the  linear  weighting  of  SSTMC  and  Tfieid  is  also  validated  by  the  dashed  curves  on  the  rightmost 
plots  in  Figures  10  through  12,  which  represent  the  values  of  100  cfield  as  a  function  of  ITfield  -  Tbuoyl.  Recall  that  the 

linear  weighting  expressed  in  Equations  (12)  and  (13)  resulted  in  the  relationship  between  the  coefficient  of  Tfield  and  the 
rms  errors  of  SSTMC  and  Tfield: 


Tfield  _ 
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G 


2 


2  2 
° Tfield  +  ^SSTMt 


(15) 


In  addition,  values  of  1.0,  0.7  and  0.55  were  the  rms  errors,  (Jrfield ,  quoted  earlier  for  Clim,  K100  and  K10,  respectively 
[2].  Values  of  crSST  will  be  visually  estimated  from  Figures  10  through  12,  as  we  now  consider  the  least  aggressive 
maximum  value  of  5.0  for  ITfield  -  Tbuoyl. 

From  Figure  10,  the  SSTMC  rms  error  is  visually  estimated  to  be  0.55.  Combining  with  the  three  characterizations  of 
Tfieid,  we  arrive  at  values  for  100  cfield  of  23.22%,  38.2%  and  50.0%  for  Clim,  K100  and  K10,  all  representing  a 
difference  of  less  5%  compared  to  values  from  the  plot.  From  Figure  11,  the  SSTMC  rms  error  is  close  to  0.75  when 
ITfield  -  Tbuoyl  <  5°K.  Combining  with  the  three  characterizations  of  Tfieid,  we  arrive  at  values  for  100  c^ield  of  36.0%, 

53.4%  and  65.0%  for  Clim,  K100  and  K10,  all  representing  a  difference  of  less  15%  compared  to  values  from  the  plot. 
Finally,  from  Figure  12,  the  SSTMC  rms  error  is  estimated  to  be  0.67  when  ITfieid  -  Tbuoyl  <  5°K.  Combining  with  the  three 
characterizations  of  Tfieid,  we  arrive  at  values  for  100 cfeld  of  30.9%,  47.8%  and  59.7%  for  Clim,  K100  and  K10,  all 
representing  a  difference  of  less  10%  compared  to  values  from  the  plot. 

The  gradual  trend  for  100  cfield  for  all  zenith  bands  is  to  increase  to  100%  as  the  aggressiveness  of  pre-filtering  increases. 
The  smoothness  of  the  dashed  curves  are  additional  evidence  of  the  inherent  stability  of  the  TfiedSST  algorithm.  This 
trend  to  a  value  of  100%  is  revealed  by  the  equation  for  cfeld  above  by  setting  the  Tfield  rms  error  to  0.0,  a  condition 

which  results  from  the  artificially  induced  reduction  in  rms  error  due  to  more  aggressive  pre -filtering.  Discrepancies 
between  calculated  values  of  100  cfield  and  those  read  from  the  plots  in  Figures  10  through  12,  could  have  arisen  from  the 

fact  that  all  the  runs  were  performed  by  allowing  all  of  the  MCSST  coefficients  to  interplay  freely  with  the  Tfieid 
coefficient.  The  MCSST  algorithm  was  not  inserted  as  a  single  predictor  by  itself  during  our  investigation,  thereby 
slightly  diminishing  the  dichotomy  between  SSTMC  and  Tfieid  needed  to  derive  the  above  relationships  between  variances. 

It  is  also  probable  that  some  small  correlation  exists  between  errors  in  SSTMC  and  Tfield  with  respect  to  Tbuoy,  which 
would  slightly  affect  the  predictions  of  the  values  for  100  cfdd  •  Figure  13  shows  the  scatter  plot  of  the  errors  in  SSTMC  - 

Tbu0y  versus  errors  in  Tfieid  -Tbuoy  using  Tfield  =  K10,  with  filtering  of  ITfieid  -  Tbuoyl  <  0.7°K  over  the  full  swath  from  0°  to 
70°.  As  was  done  previously,  both  the  plotting  sequences  from  Blue-to-Red  and  Red-to-Blue  are  displayed  due  to 
overlap  between  zenith  bands.  The  rather  negligible  correlation  value  of  less  than  0.20  supports  our  earlier  assumption  of 
uncorrelated  errors  that  resulted  in  Equation  (16),  an  important  relationship  which  guarantees  that  the  rms  error  of 
SSTTfield  will  always  be  less  than  the  rms  error  with  respect  to  Tbuoy  of  either  SSTMC  or  Tfieid. 

Overall,  the  interpretation  of  the  TfieldSST  algorithm  as  a  simple  linear  combination  of  daily  satellite -derived  data 
through  SSTmc  and  a  priori  satellite-derived  background  through  Tfield  resulted  in  predictions  for  the  values  of  their 
coefficients  that  closely  mimic  both  the  trends  and  magnitude  of  the  results  of  our  investigations. 
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Figure  13.  Scatter  plot  of  errors  in  SSTMC  -Tbuoy  versus  errors  in  Tfield  -Tbuoy  using  Tfield  =  K10. 

We  finalize  this  discussion  with  a  simple  analogy  that  parallels  our  findings  and  neatly  summarizes  our  results,  which 
are  very  similar  to  the  well-known  noise  reduction  achieved  by  frame-averaging  imagery  captured  by  a  camera.  If  the 
noise  is  assumed  to  be  uncorrelated  from  frame  to  frame  and  of  equal  variance,  the  addition  of  only  two  frames  reduces 
the  variance  by  half  or,  equivalently,  reduces  the  noise  rms  error  by  30%,  a  substantial  result  that  was  shown  earlier 
during  our  discussion  of  linear  weighting  between  SSTMC  and  Tfield. 


If  this  “imagery”  now  represents  the  global  distribution  of  Tbuoy,  which  is  the  information  that  we  are  trying  to  extract, 
then  SSTmc  and  Tfield  are  equivalent  to  two  “noisy”  frames  of  Tbuoy,  but  with  different  variances.  As  we  found,  the  result 
of  adding  them  together  with  linear  weighting  is  to  automatically  reduce  the  variance  of  SSTTfield.  The  TfieldSST 
algorithm  finds  the  optimal  combination  of  these  two  noisy  frames  that  will  provide  the  smallest  resulting  variance, 
thereby  providing  a  less  noisy  image  of  Tbuoy. 


4.  CONCLUSIONS 

Results  from  an  in-depth  investigation  about  the  potential  benefits  of  using  Tfield,  a  first  guess  temperature  field  currently 
used  in  NLSST,  as  a  separate  predictor  in  the  MCSST  algorithm  were  presented  and  scrutinized  for  clues  that  eventually 
led  to  a  simple  linear  weighting  formulation  named  TfieldSST.  A  well-established  relationship  of  the  “resistors  in 
parallel”  form,  Equation  (16),  between  the  linear  regression  coefficients  of  the  two  predictors  of  TfieldSST,  SSTMC  and 
Tfield,  and  their  individual  rms  errors  guarantees  that  the  rms  error  of  SSTTfield  will  always  be  less  than  the  rms  error  of 
either  component.  Unlike  other  gap-free  background  temperature  fields  which  use  in  situ  measurements,  the  first  guess 
temperature  field  developed  at  NAVOCEANO  is  derived  from  previous  satellite  data,  making  the  new  TfieldSST 
algorithm  entirely  dependent  on  satellite  observations  only,  both  past  and  present. 

It  was  shown  why  a  genuine  minimum  30%  reduction  in  MCSST  rms  error  is  to  be  expected  from  the  TfieldSST 
algorithm,  thereby  automatically  increasing  its  accuracy  from  a  current  0.5 °K  to  0.35°K.  In  addition,  we  showed  how  the 
ultimate  goal  of  0.3 °K  rms  error  could  easily  be  achieved  with  the  TfieldSST  algorithm  at  the  very  minimal  cost  of  a  5% 
reduction  in  the  80  %  fraction  of  the  original  buoy  data  points  over  the  entire  range  of  zenith  angles  from  0°  to  70°, 
achieved  with  a  more  aggressive  filtering  with  respect  to  Tfield,  ITfield  -  Tbuoyl  <  0.7°K  .  This  is  to  be  compared  to  the 
present  30%  reduction  in  the  number  of  buoy  data  points  over  the  restricted  zenith  angle  range  form  0°  to  53°,  with  a 
rms  error  of  about  0.5 °K  achieved  with  the  filtering  threshold  of  ITfieid  -  Tbuoyl  <  2.0°K.  It  should  be  emphasized  that 
although  Tfield  is  presently  used  in  the  NLSST  algorithm,  its  information  content  is  “tamed”  due  to  the  fact  that  it  comes 
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in  as  a  multiplier  to  another  predictor,  Tn-Ti2.  In  sharp  contrast,  the  TfieldSST  algorithm  makes  the  optimal  use  of  Tfield 
information  by  using  it  as  a  separate  predictor. 

The  simple  linear  weighting  form  of  the  TfieldSST  algorithm  inevitably  puts  Tfieid  on  equal  footing  with  daily  satellite 
retrievals,  so  that  minimization  of  rms  error  is  now  two-pronged:  improvements  to  SST  algorithms  and  improvements  to 
the  accuracy  of  Tfieid.  Efforts  are  underway  at  NAVOCEANO  to  produce  a  K2  first  guess  temperature  field  that  should 
result  in  reduction  in  its  rms  error.  Increasing  the  accuracy  of  existing  SST  algorithms  at  large  zenith  angles  beyond  53° 
should  be  achievable  by  including  sea  surface  emissivity  and  roughness,  as  well  as  hardware -centric  reduction  of  rms 
error  with  the  acquisition  of  multiple  frames  of  brightness  temperatures  for  the  same  “pixel”. 

The  findings  and  conclusions  of  this  study  were  limited  to  a  NAVOCEANO -generated  buoy  data  set  and  VIIRS  satellite 
data  from  June  2012  and  will  therefore  require  examination  of  seasonal  trends  for  algorithm  stability.  However,  we  can 
conclude  that  the  generic  form  of  the  TfieldSST  algorithm  predicts  that  the  addition  of  Tfield  to  any  existing  SST 
algorithm,  both  daytime  and  nighttime,  will  result  in  significant  reductions  in  rms  error  over  the  full  swath  of  zenith 
angles  and  enhance  present  operational  capability  for  global  coverage.  We  therefore  strongly  recommend  that  Tfield  be 
included  in  the  VIIRS  EDR  data  stream  for  use  by  developers  of  SST  algorithms. 
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